Temperature-dependent terahertz conductivity of topological insulator 
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Using Terahertz Time-Domain Spectroscopy, we study the temperature-dependent complex opti- 
cal conductivity of the topological insulator, Bi1.5Sbo.5Te1.sSe1. 2 single-crystal from 5 K to 150 K 
in the terahertz regime (0.4 - 3.0 THz). We analyze our experimental results using the Drude- 
Lorentz model, with the Drude component representing the metallic surface state and the Lorentz 
term representing the bulk insulating state. We find the conductivity to be dominated by the Drude 
contribution, suggesting the presence of metallic surface states. The low-frequency real conductivity 
follows a thermally-activated behavior. Its origin is also discussed. 



Topological insulators are electronic materials with 
an insulating bulk and conducting surface [l|, [2|. Such 
unique features have been attributed to the presence of 
time-reversal symmetry and spin-orbit interaction 0, Q . 
Bi2Se3, Bi2Te3 and Sb2Te3 are several of the semi- 
metallic chalcogenides shown to have the properties of 
a three-dimensional topological insulator [4H6fl. Inves- 
tigations of the surface state have been made by sur- 
face probes such as Angle- Resolved Photoemission Spec- 
troscopy (ARPES) H 0, [|| and Scanning Tunneling 
Microscopy (STM) 0, Q. In Bi 2 Se 3 and Bi 2 Te 3 , the 
presence of non-stiochiometry-induced bulk carriers in- 
duced extrinsic metallic conduction in the bulk (Tol — 
13|, making it difficult to detect the metallic surface 
states. Recently, new topological insulators Bi2Te2Se 
and Bi1.5Sbo.5Te1.7Se1. 3 were found [l4, 15 1, where non- 
stiochiometry-induced donors and acceptors compensate 
each other, yielding high bulk resistivities of 1-10 il.cm 
and high contribution (~70%) of the total conductance 
by surface transport, thus showing a greater contrast 
between the bulk and surface resistivities compared to 
Bi 2 Se 3 , Bi 2 Te 3 and even Bi 2 Te 2 Se ARPES exper- 

iments have also shown that their chemical potential is 
always located inside the bulk band gap (l7l ). 

Terahertz Time-Domain Spectroscopy (THz-TDS) is 
a non-contact far-infrared optical technique suitable for 
probing the low-energy excitations of strongly correlated 
electronic systems such as cuprate superconductors [18j . 
pnictide superconductors [ijj], and colossal magnetoresis 



tance manganites [18] . Recently, Aguilar et al. |20l. |21| 
applied the same technique to study the THz response 
of topological insulator Bi2Se3 thin films, and had to ac- 
count for free carriers from the bulk in characterizing the 
conductivity of the material. In this Paper, we present 
THz-TDS studies of Bii. 5 Sbo. 5 Tei. 8 Sei.2 (BSTS) single 
crystal at temperatures 5 K-150 K, to study the fre- 
quency and temperature-dependent optical conductivity 
in the far-infrared regime (0.4 - 3.0 THz). By model- 
ing the total conductance as the sum of Drude (surface) 



and Lorentz (bulk) components, we obtained the surface 
conductivity <r sur f ~10 5 (fi.cm) -1 to be 10 4 times the 
bulk conductivity <7buik ^lO 1 (Sl.cm) -1 . This suggests a 
significant contribution of the topological surface states 
to the conductivity of the BSTS sample. 

Bi1.5Sbo.5Te1.8Se1. 2 single crystals were synthesized 
by melting high-purity (99.9999%) Bi, Sb, Te and Sc 
with molar ratio 1.5:0.5:1.8:1.2 at 850°C in an evacuated 
quartz tube. The temperature was then gradually de- 
creased to room temperature over a span of three weeks 
[2^ | . The BSTS single crystal was then cleaved repeat- 
edly to a thickness of ^60 /im while maintaining a large 
surface area (~4x4 mm 2 ). 

THz transmission of the BSTS single crystal were 
measured using a conventional THz-TDS system (Ter- 
aView Spectra 3000) incorporated with a Janis ST- 100- 
FTIR cryostat. The THz signal was generated and de- 
tected by photoconductive antennae fabricated on low- 
temperature-grown GaAs films. The aperture diameter 
is 3.5 mm, allowing for an accurate measurement of the 
THz signal down to ^0.4 THz. The time-domain electric 
field of the THz pulse signal is transmitted through the 
BSTS sample (Es(t)), while the reference signal (E^(t)) 
is transmitted through vacuum. 1800 THz traces were 
taken in 60 seconds for each reference or sample run. 
The sample holder was moved back and forth between 
the sample and reference positions by means of a verti- 
cal motorized stage with a resolution of 2.5 ^m. Fast 
Fourier Transform (FFT) was then performed on the 
time-domain THz signal to obtain the amplitude and 
phase of the THz spectra. Since the THz-TDS detects 
both the amplitude and phase of the THz signal, there is 
no need to use the Kramers-Kronig transformation to ex- 
tract the real and imaginary components of the material 
optical parameters. 

By finding the ratio between the sample [£^5(0;)] and 
reference [£^(0;)] spectra, one is able to obtain the 
complex transmittance, T{uj)=Es(w)/ Er{lS). Figure [1] 
shows the amplitude of the experimental transmittance 
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FIG. 1: Amplitude of complex transmittance, |T(cj)|, at var- 
ious temperatures. 



of BSTS, |T(u;)|, at various temperatures. A series of 
wiggles have been observed at low temperatures which 
eventually get washed out with increasing temperature 
above ~100 K. We associate these wiggles to the pres- 
ence of impurity states in our BSTS sample due to dop- 
ing. Similar wiggles at low temperatures were also seen 
in optical conductivity data of Bi2Se3, Bi2Se2Te, which 
the authors associated with hydrogen-like Is— >np tran- 
sitions of isolated impurities [16]. In addition, a drop in 
\T(ui)\ to zero has been observed at ~1.7 THz and af- 
ter which, rises again at ~2.8 THz. This corresponds 
to an optical phonon mode at 1.9 THz, to be described 
later. By fitting the experimental T(u>) to the theoretical 
expression 
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one is able to able to obtain the complex refractive index 
of the sample, n = n+in. Here d (= 60 /im) is the sample 
thickness, c is the speed of light in vacuum, and M=6 
is the number of multiple reflections of the THz pulse 
within the sample. Note that the wiggles seen below 80 K 
are not a result of multiple reflections within the sample 
— their existence is robust when we vary M around its 
expected value. 

The complex refractive index is then used to calcu- 
late the complex conductivity a = o\ + i<72, where 
<7i(w) = 2nKLueo and a^iS) — (fee — n 2 + k 2 )lo£q, eo being 
the permittivity of free space, and the high-frequency 
dielectric constant. As an unknown quantity, was 
initially set to 1 when we calculate the experimental con- 
ductance, via the expression Gdata{u) = &(ui)tbuik', while 
it will be used as a fitting parameter [20] in the conduc- 
tance fitting for each temperature. Since the bulk state 
is highly insulating compared to the conducting surface 
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FIG. 2: (a) Real conductance Gi(cj) at 30 K. Inset shows 
the imaginary conductance GiifJ). (O) = data. Solid lines = 
Drude (green) and Lorentz (blue) contributions, and the fit 
to data (red), (b) Low- frequency real conductivity a\ versus 
temperature. (O) = data. Solid line = fit to Eq. ([3|. 



states, the transport dynamics of the surface states may 
be modeled by a Drude (D) term, whilst the insulating 
bulk state by a Lorentz (L) term. Therefore, we model 
the conductance as G mode i(u) = a mo del(^)tbulk, where 

€qUJ 2 jj 

o mo dei(u) =a D (u})f + a L (uj) = — / 
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with / = t sur f /tfruik being the contribution of the 
Drude conductivity to the model conductivity. The sur- 
face state thickness t sur f^2 nm is taken to be that 



of high carrier-density Bi2Se3 [23j], and bulk thickness 
tbuik=d>0 fim. The fitting parameters lo p d and jd de- 
note the plasma frequency and scattering rate of the 
Drude component, respectively. wol and jl are 

the plasma frequency, oscillator frequency and the scat- 
tering rate of the Lorentz component, respectively. The 
temperature dependence of these fitting parameters will 
be discussed later. Figure [2ja) shows the real and imag- 
inary components of the experimental complex conduc- 
tance G(uj) = Gi(uj) + zG*2(w) of the sample at 30 K, 
alongside the fit to the Drude-Lorentz model, and the 
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respective Drude and Lorentz contributions. Notice that 
the dips in transmission at 0.7 THz, 0.9 THz and 1.2 THz 
in Fig.Q]become peaks in G\ and at the same frequen- 
cies, in Fig. HJa). 

We also tried to fit our data using two Drude and one 
Lorentz terms, where the extra Drude term representing 
the free carriers in the bulk, as has been done in a THz- 



TDS study of Bi 2 Se 3 thin films [2l|. We found that the 
Drude fitting parameters are now not well-constrained 
and depend sensitively on initial guesses. In Bi2Se3 
thin films, a strong Drude component was observed at 
sub-THz frequencies and requires both bulk and surface 
Drude terms to account for it [2l[. The low- frequency 
response of our BSTS sample, on the other hand, is 
Lorentz- like with a Drude offset (see Fig. Ufa)). This 
is consistent with our BSTS sample having a large bulk 
resistivity, thus needing only one Drude term (represent- 
ing the surface) in our conductance fit. 

We estimate how dc conductivity of BSTS varies 
with temperature by plotting the low-frequency 
(Ti (w=0.4 THz) versus temperature, then fitting it to 
the thermally-activated hopping model given by [24j 



o-(T) = Aexp [-A/fcuT] + D 



(3) 



from 5 K to 150 K, where A is a constant, D is a non- 
zero offset at K, A is the activation energy, and ks 
is the Boltzmann constant. Figure EJb) shows the low- 
frequency conductivity together with the fit to Eq. ([3]), 
with D — (2.80 ± 0.01) (O.cm) -1 , and activation energy 
A = (26.7±0.2) meV. The fitted value of A is consistent 
with the results by Ren et al. that vary between 22 and 
53 meV [2^. A similar energy scale of 30-40 meV has 
also been found in Bi2Te2Se, which has been attributed 
to transitions from the impurity bound states to the elec- 
tronic continuum (l6| . 

Figures [3](a) and (b) show the real Drude and Lorentz 
conductance (Gf rude and G\ orentz ) and conductivity 
(of rudc and crL° rentz ) at selected frequencies. cf rude and 
^Lorentz nave been obtained by dividing each conduc- 
tance component with the surface and bulk thickness 
respectively. At a fixed frequency, the Drude conduc- 
tance increases with increasing temperature, but at a 
fixed temperature, it decreases slightly with increasing 
frequency. On the other hand, the Lorentz conductance 
is temperature-independent, but increases significantly 
with increasing frequency. The frequency dependences 
of the Drude and Lorentz components are just natural 
consequences of the respective models at low frequencies. 
Figure shows that the surface (Drude) conductivity 
Csurf ~10 5 (fLcm)" 1 is ~T0 4 times the bulk (Lorentz) 
conductivity Ubuik ~10 (fi.cm) -1 . This suggests a sig- 
nificant contribution of the topological surface states to 
the conductivity of the BSTS sample. Note that the 
conductance values are independent of t sur f and t^uiki 
whereas the surface (bulk) conductivity value scales in- 
versely as t sur f (tbuik). If we take t sur f to be 40 nm as in 
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FIG. 3: Temperature dependence of the (a) Drude and (b) 
Lorentz, conductance and conductivity, at selected frequen- 
cies. 
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FIG. 4: Temperature dependence of the (a) Drude spectral 
weight ujp D , (b) Lorentz spectral weight ujp L , (c) Drude scat- 
tering rate 70, and (d) Lorentz scattering rate 7_l. Solid line 
in (a) is fit to Eq. ©. 



low carrier-density Bi2Se3 [13] . then the surface-to-bulk 
conductivity ratio will be one order of magnitude less, 
which is in better agreement with electrical transport 
measurements |14| . Furthermore, if we were to replace 
tsurf with tbuik, the resulting surface-to-bulk conductiv- 
ity ratio (~10 _1 -10°) would have been unphysical. 
We next examine the factors that contribute to the fre- 
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quency and temperature dependence of the surface and 
bulk conductance, by looking at the temperature depen- 
dence of the fitting parameters in Eq. Figures Hta) 
and (b) display the temperature-dependent Drude spec- 
tral weight (ujp D ) and Lorentz spectral weight (oJp L ) re- 
spectively, while Figs.rjjc) and (d) show the surface (jd) 
and bulk (jl) scattering rates respectively. We found 
uip D to increase significantly with temperature, while 
7d (proportional to charge carrier mobility) shows little 
change with temperature. As for bulk carrier dynamics, 
ujp L exhibits insignificant correlation with temperature 
while 7£ shows a slight increase with temperature. Since 
spectral weight is directly proportional to the charge car- 
rier density N via the relation uj 2 — Ne 2 /eom*, where e 
denotes the elementary charge and m* is the effective 
mass of the charge carrier, this suggests the increase in 
the surface charge density with temperature while the 
scattering rate changes only marginally. This is unlike 
the case of conventional metals, where the increase in re- 
sistivity with increasing temperature originates from the 
increase in scattering rate, with the charge density re- 
maining temperature independent. The increase in sur- 
face charge carrier density may be due to the contribu- 
tion of bulk carriers from the donor impurity band to 
the topological surface states (Fig. [5]). This band, sit- 
uated ~30 meV (=activation gap A) below the Fermi 
level, starts to contribute carriers to the Dirac cone above 
60 K, leading to a rise in the Drude spectral weight uj 2 d . 
This relationship between A and u> 2 D is supported by 
our fit of uJ 2 pD {T) (in Fig. Ha)) with Eq. © — a good 
fit was obtained with A=(31 ± 4) meV (solid line in 
Fig. IHa)), consistent with the earlier-obtained value of 
A=(26.7±0.2) meV from the ct(w=0.4 THz) fit. This 
is not surprising, since, from Eq. ©, G(co 0,T) oc 
pX)(T). Note that the temperature dependence of oj 2 d , 



and hence value of the activation gap A, do not depend 
on the value of t sur f . Figure [5] shows a schematic of the 
relative positions of the valence band, Dirac point, im- 
purity band, chemical potential, and conduction band in 
BSTS. From oj 2 d , one obtains the carrier density of the 
Drude (surface) term N sur f~l0 21 cm~ 3 . Compare this 
with the carrier density of a typical metal of ~10 23 cm -3 

From the conductance fittings, we also obtained the 
oscillator frequency wol = (1.92 ± 0.04) THz, and tem- 
perature independent. This oscillator can be attributed 
to an optical phonon mode, which also agrees with the 
THz-TDS study by Aguilar et al. on Bi2Se3 thin films 
[2ll |. where a similar optical phonon has been found at 
~2.0 THz. This result is also consistent with the A] lon- 



gitudinal optical phonon obtained by pump-probe 27I.|28| 
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FIG. 5: Schematic of band diagram of BSTS, with the 
(donor) impurity band lying just below the Fermi level. The 
relative positions of the valence band, Dirac point, Fermi 
level and conduction band were taken from Ref. [FT] for a 
Bi1.5Sbo.5Te1.7Se1. 3 sample. A= activation gap. 



of Bi1.5Sbo.5Te1.8Se1. 2 as a function of temperature using 
THz-TDS. By modeling the experimental conductance as 
a sum of a Drude term (from the surface) and a Lorentz 
term (from the bulk), we found the surface conductivity 
to be 10 4 times greater than the bulk conductivity. The 
high surface conductivity compared to the bulk suggests 
the presence of metallic surface states of BSTS. 
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In conclusion, we measured the terahertz conductivity 
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